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Abstract 


Visualisation of the optical field characteristics within the source-fiber region in transmitter 
packages is presented using the Finite DiiFerence-Time Domain (FD*-TD) method in two 
dimensional configuration. The reflection characteristics of the optical field from the fiber 
end face and the coupling efficiency are studied for three coupling arrangement schemes viz. 
centering a flat fiber end face directly over the light source; using a taper-end fiber; and 
placing a naicrosphere in between the light source and the flat fiber end face. Finally, the 
effect on the field pattern and the coupling efficiency due to the misalignment is studied only 
for the third coupling scheme. 


1 



Chapter 1 


Introduction 


1.1 General 


The main sections of an optical fiber communication link (Figure 1.1) are a transmitter 
consisting of a LED or laser diode and its associated drive circuitry, an optical fiber as the 
transmission medium, and a receiver consisting of a photodetector plus an amplifier and 

TRANSMITTER 

j 

I 

: 



Figure 1.1: Basic configuration of an optical fiber transmission link. 
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signal- restoring circuitry. In the transmitter optical power is coupled from the source into 
the fiber. 

When the source is turned on optical field propagates towards the fiber and a part of the 
input energy is reflected back towards the source at the air-core and air-cladding boundaries, 
and a part of the remaining portion of the input energy couples into the fiber and propagates. 
At the receiving end of fiber, power propagates towards a photodetector in the receiver which 
converts the optical energy into the electrical signal. 

A measure of the amount of optical power emitted from a source that can be coupled into 
a fiber is given by the coupling efficiency tj t=-^ where Pf is the power coupled into the 
fiber and Pg is the power emitted by the source. The more the coupling efficiency the better 
the performance of the communication link. The coupling efficiency depends on the optical 
source used, the type of fiber and on the couphng process, whether or not lenses or other 
coupling attachments are used. 

Nowadays, people are interested in the packaging of optical systems. For a compact 
packaging of the optical transmitter/receiver systems, we must about the field characteristics 
between the optical source/ detector and the coupling end of the fiber. It is also important 
to know about the coupling efficiency for different attaachment schemes. Visualisation of the 
optical field characteristics within the narrow coupling regions mentioned above is neither 
possible experimentally nor theoretically. The main purpose of this thesis is to visualize the 
optical fields in the coupling regions of transmitter packages and to throughly determine the 
improvements, if any, due to coupling attachments and degradations in coupling owing to 
misallingments. 

1-2 Some Coupling Models 

Several possible coupling configurations are given in Figures 1.2 to 1.9 [16, 17, 18, 19, 20, 21, 
22, 23, 24, 25]. In Figure 1.2 we have shown optical power-launqbing based on centering a flat 
fiber end face directly over the light source and it is kept as close as possible to the source [29]. 
In this configuration, if the source emitting area is larger than the fiber core area, then 
the resulting optical power coupled into the fiber is the maximum that can be achieved. In 
Figure 1.3 a taper-ended fiber is centered instead of the flat fiber end face [21]. In Figure 1.4 
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Source 

Figure 1.2: Power coupling based on centering a flat fiber end face directly over the source. 




Figure 1.3: A taper-ended fiber is centered over the light source. 

Nonimaging 

microspliere 




Figure 1.4: A noniinaging microsphere is placed between source aaid fiber end face. 
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Active 



Source 



Figure 1.5: A larger spherical lens is used to image the source on the core area of the fiber 
end. 



Figure 1.6: Coupling configuration using gradient-index lens. 


Cladding 



Figure 1.7: A rounded-end fiber is used in place of a flat fiber 
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Cylindrical lens 




Figure 1.8: A cylindrical lens formed from a section of jfiber is used in between tbe source and 
the flat fiber end. 



Figure 1.9: Coupling configuration consisting of a spherical-surfaced LED and spherical-ended 
fiber. 





a nonimaging microsphere [23] is placed in between the source and the flat fiber end face. 
This is kept in contact with both the fiber and the source. In Figure 1.5 a larger spherical 
lens is used to image the source on the core area of the fiber end [20]. In Figure 1.6 we have 
shown the coupling configuration using graded-index lens [16, 17, 18, 19, 20]. In Figure 1.7 
a rounded-end fiber is used in place of a flat fiber end face [22]. A cylindrical lens generally 
formed from a short section of fiber is placed in between the source and the flat fiber end face 
as shown in Figure 1.8 [29]. Finally, in Figure 1.9 we have presented a coupling configuration 
which consists of a spherical- surfaced LED and spherical-ended fiber [25]. 

Although these techniques can improve the coupling efiiciency over the butt coupling (Fig- 
ure 1.2), they also cause additional fabrication and handling difficulties. In case of taper-ended 
fiber the mechanical alignment must be done with greater precision since the coupling effi- 
ciency becomes a more sharply peaked function of the spatied alignment. However, alignment 
tolerances are increased for other types of lensing schemes [21]. 

Now some questions arise: what are the optical field characteristics within the regions 
between the source and the fiber end and how can we visualize these field characteristics for 
different configurations of coupling? How does the coupling efficiency change for diflferent 
configurations of coupling as shown in Figures (1.2- 1.9)? How does the coupling efficiency 
vary due to the misalignment for the above configurations? 

The first configuration (Figure 1.2) is the simplest one. The second (Figure 1.3)and third 
(Figure 1.4) configuratios are somehow more compliacted than the first configuration. These 
three configurations are very easy to implement in the computer. That is why, we consider 
the above three as the basic configurations to answer the above questions. 

1.3 Problem Description 

If we take a laser as the light source, its output intensity distribution in the transverse 
directions is Gaussian when it oscillates in the fundamental mode. The output of the source 
can be given by a set of electromagnetic field components. When the source is turned on 
electromagnetic waves propagate towards the fiber coupling end and are diffracted. But, 
the diffraction formulas available in the literature are the approximations and these are only 
applicable beyond a certain distance from the source. Therefore visualisation of the optical 
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field characteristics is not possible analytically within a small distance(in jj-m range). That is 
why we need numerical methods to visualize the above field characteristics. 

In this thesis we have tried to visualize the optical field characteristics within the regions 
between the Source-fiber in the transmitter system using a computational technique named 
Finite Difference Time Domain (FD-TD) method [1, 2, 3, 4]. We have also tried to com- 
pare the coupling efficiency for three coupling models (Figure ??- ??). Though there are 
several methods like TLM (Transmission Line Matrix) [15], Moment method [14], Iteration 
method [12] etc. to solve the field problems numerically, we have used FD-TD method due to 
its high accuracy and ease of implementation on computers. By this method the Maxwell’s 
equations are discretized in both space and time and then solved in time marching sequence 
by alternatively calculating the electric and magnetic fields in an interlaced grid. Boundary 
condition named absorbing boundary condition (ABC) should be applied on the boundaries 
of the computation domain to simulate the unbounded physical surroundings. We have used 
Mur’s second order ABC to implement the FD-TD method [5]. 

To implement the FD-TD method in three dimensional configuration, it requires a large 
amount of computer memory as well as computation time. Therefore we have studied in the 
two dimensional configurations only. We have assumed that the fields do not vary in the x 
direction. Fiber is assumed as the planner slab wave guide with axis along the z axis. We 
have studied only for TE mode of propagation with field components E^, Sy and Hz- Other 
field components are assumed to be zero. Our study can be easily extended to TM modes. 

As the equilibrium numerical aperture of a fiber can be obtained after propagating a long 
distance which may be several meters (for long fibers it is taken as 50 m), an equilibrium 
mode can not be obtained by analyzing the FD— TD method. Therefore, the study of the 
coupling efficiency accurately is not possible. 

We have studied only three coupling configurations given by Figures 1.2 to 1.4 for single 
moded step index fiber with core width and cladding width=120^m(total),‘ core 

refractive index=1.50 and cladding refractive index= 1.499. 
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1.4 Overview Of The Thesis: 


In. chapter-2 we present a brief description of the Finite Difference-Time Domain method 
and Mur’s absorbing boundary condition (ABC). In chapter-3 the results and discussions are 
given for three transmitter packages. The further scope of work is discussed in chapter-4. 
Finally, the conclusions are presented in chapter-5. 
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Chapter 2 


The Finite Difference Time 
Domain Method 

2.1 Introduction 

In this chapter, we give a brief description about the Finite Difference - Time Domain (FD 
- TD) method [1, 2, 3, 4] and Mur’s absorbing boundary condition (ABC) [5] which has 
been used to solve our thesis problem. The FD - TD method is a very useful method to 
solve Maxwell’s curl equations numerically and it is a very handy computational technique 
for a wide variety of problems in electromagnetic radiation, scattterering and diffraction in 
both homogeneous and inhomogeneous media. Simulation of electromagnetic problems using 
the FD - TD method was first proposed by K.S.Yee in 1966 [1]. Following Yee’s method 
Maxwell’s equations are discretized both in time and space and then solved in time marching 
sequence by alternatively calculating the electric and magnetic fields in an interlaced grid. 
The resulting FD - TD method allows propagation of waves in all directions. 

Many of the applications involve modeling electromagnetic fields in an unbounded space. 
But, due to the limited storage space of the computers, numerical computation space must be 
finite. Therefore, a certain kind of boundary condition, which is called the absorbing boundary 
condition (ABC), needs to be applied on the artificial boundaries of the computation domain 
to simulate the unbounded physical surroundings. Though there are several higher order 
ABC developed so far in the literature [6, 7, 8, 9] Mur’s second order ABC is used all through 
this thesis to simulate the infinite physical space to a finite computation domain due to its 
less amount storage requirements and ease of implementation. 

In section 2.2, a brief review of the derivation of the discretized form of the MaxweE’s 
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equations is presented. In section 2.3, Mur’s ABC is given in discretized form. In section 2.4, 
software implementation of the FD— TD method to study the propagation characteristics of a 
point source and some numerical results using the above ABC are presented. Finallly, section 
2.5, gives some concluding remmarks and discussions. 

2.2 Maxwell’s Equations And The Equivalent Set Of Finite 
Diifference Equations 


The relations among the electric field vector (E), magnetic field vector (H), electric displace- 
ment vector (D), magnetic flux density vector (B) and current density vector (J) are given 
by the Maxwell’s equations; 




0 
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where /x, e, cr, p are the permeability, permittivity, conductivity and charge density respec- 
tively and are assumed to be function of space and time. Assuming the medium to be charge 
free (i.e. p = 0) and using rectangular coordinate system, (2.1) and (2.2) can be expressed to 
the following set of scalar equations: 
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Figure 2.1: Positions of th.e field components about a unit cell of the'Yee lattice. The E - 
components are in the middle of the edges and the H - components are in the center of the 
faces. 

where the subscripts x, y and z denote the x, y and z components of E, H, D, B, and J 
vectors and each component is a function of x, y, z and t. 

We denote a grid point of the space as 


{i,j,k) = (iSx,jSy,kSz) (2.14) 

and for any function of space and time we put 


F(iSx,j6y,kSz,n6t) = F^(i,j,k) (2.15) 

Now, if we use a centered difference scheme for both space and time, eqn. (2.8) can be dis- 
cretized to : 


-f |,fc -f |) - J9a; + |;fe+ I) 

St 

_ + |, A; -H) - E^{i,j + \,k) 

Sz 

E^(i,j + l,k+ I) - E^{i,j,k+ 1) 

Sy 

For (2.11) we have 


(2.16) 
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- DT^{i + 

6t 

_ Hz +h^) - 

Sy 

nrhi + l.i, k) - nrhi + hJ, k - 1) 

Sz 

+ jrHi+\j,k) 


(2.17) 


Tlie above discretization rule to solve the electromagnetic field comonents numerically, was 
first proposed by K.S.Yee and therefore it is also called the Yee’s algorithm. According to 
this algorithm, the positions of the electric and magnetic field components can be shown as 
in Figure 2.1 . 

Using the same discretization rule for other scaler equations in (2.8)-(2.13) and (2.3)- 
(2.7); rearranging the terms and considering the programming constraints, the electric and 
magnetic field components can be given in discretized form by the following six equations [11]: 


= 1 - 
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( 2 . 21 ) 
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( 2 . 22 ) 
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Prom tie above algorithm we see that, if the components of E are calculated at nSt, where 
St is the time step, and n is any nonnegative integer, the components of H are calculated at 
(n + \)St. For this reason, this algorithm is also called the leapfrog method. 

For any finite difference scheme, a stability condition must be found which guarantees 
that the numerical error generated in one step of the calculation does not accumulate and 
grow. To ensure the stability of the above time - stepping algorithm, St is to be chosen to 
satisfy the following condition: 


'^max — 





(2.24) 


where Vmax is the maximum wave phase velocity expected within the medium. 


2.3 Mur’s ABC 


Mur’s second order ABC for a; = 0 boundary plane can be given in the form of the following 
differential equation; 


'IJL. 

Vp dx dt 


l/_^ 



(2.25) 


where U can be any electric field component at the boundary x — 0 and Vp is the phase 


velocity of the wave in the medium. A popular discretization of 2.25 is 

Vp ^fJp 

+i [didI. {vSjj. + + OiDL ■ = 


(2.26) 


Here, Di and are the standard forward, backward and central difference operators 
and Uf'j k = U(iSxJSy, kSz, nSt); Sx, Sy and Sz are the grid sizes in the x, y and « directions, 
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respectively, and 6t is the time step. The above equation (2.26) can also be given as 




Tjn~\ , — 8 

Vp6t + 8 


{vm + psji) 
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(2.27) 


where we have assumed that 8x= 8y = Sz = 8, The boundary conditions for other planes of 
the artificial boundaries can be derived in the similar way as mentioned above. 

For two-dimensional electromagnetic - field problems, the ABC of the second order ap- 
proximation can be given in a very simple form. Suppose, the fields do not depend on z and 
are E - polarized, i.e.,E = Eziz and H = Hj-ix +Hyiy . Now, (2.25) applies to Ez only and 
in this case we have 




8H. 


8Ez 


(2.28) 


8t 8y 

Substituting (2.28) in (2.25), putting 8Ezl8z = 0, with U =Ez, integrating with respect to t 
and using Ez=0 for t< 0 we get 

f 8Ez 1 8Ez jiVp 8Hx\ 


8x Vp 8t 


2 8y 


= 0 


07=0 


The discretized form of (2.29) can be given as follows: 

. - 

flVpSt frr'^i-kfn - i i\ rr^+l/'n ‘N 

^ ' Hx 2(0,J-f l)-5^s “(0,j) 




-H 5r^(i,i-i-i)-.ffr^(i-,i)) 


(2.29) 


(2.30) 


We see that (2.29) is much simpler than (2.25) but still of the same order of approximation. 
A similar ABC for the case of H - polarization can easily be derived. 


2.4 Software Implementation Of The FD— TD method: Study- 
ing With A Point Source 


la this section we give the program outline for the FD~TD method to study the field charac- 
teristics of a point source in two dimensional configuration and for an isotropic^ homogeneous, 
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source-free and nonmagnetic medium (i.e. air) in the y-z plane. The whole program can be 
given by the following steps. 

1. Begin: 

• INPUT PARAMETERS: 

2. READ: number of nodes in y-direction (ymax), number of nodes in z-direction (zmax), 
operating wavelength (A), dielectric constant of the medium (e) and permeability of the 
medium (fio)- 

• SELECTION OF SOME PARAMETERS FOR THE FD-TD METHOD: 

3. Calculate: phase velocity in the medium (vp)=^^. 

4. Select: 6y = 6z = S = -^ and St = 

5. Set the initial values of all the field components equal to zero. 

6. Select: number of time steps (umax) > ”iaa: < 2 ymax,^ ^max > 

7. Select €(j, k) = € for all values of j and k in the computation domain. 

• STARTING OF THE CALCULATION OF THE FIELD COMPONENTS: 

8. Start time stepping (n) at n = 0 

9. Choose a point (jo,ko) (say). Apply the input at this point a.s given by the equation 
(2.36) . 

• CALCULATION OF THE FIELD COMPONENTS AT THE INTERRIOR NODES: 

10. Calculate the field components at the interrior nodes of the computation domain with 
the regular finite difference scheme. 

• SPECIAL TREATMENTS OF THE BOUNDARY NODES 

11. Apply Mur’s second order boundary condition to calculate the electric field components 
at the boundary nodes of the computation domain. 

12. Increase n by one i.e. n 

13. Verify: if n < nmo® go to step ( 8); else finish. 
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Start 



Figure 2.2: Flow chart for studying with, a point source by the FD~TD method 
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14. Print outputs 


15. End 


Tlie steps are shown in the flowchart (Figure 2.2) 
2.4.1 Numerical Results 


In this section numerical simulation is performed to test the accuracy of the ABC considered 
in the previous sections. We shall present results for two - dimensional configurations only. 
We present results for fields that do not depend on x and are TE modes i.e., E = Ea-ix and 
H = Hyiy + Hziz'i where ix, iy and iz are the unit vectors along the x, y and z directions, 
respectively. A square two - dimensional mesh is used in the y-z plane with a size 76 by 76 
i.e., 0 < j < 76 and 0 < fe < 76. Absorbing boundary conditions are applied on all four sides 
of the mesh. A monochromatic isotropic point source of wavelength X=1.5iJ,m is assumed to 
be turned on at t=0 in air and 6 = Sy = Sz =2co6t = where cq is the velocity of light in 
free space. Actually,the value of Sy, Sz and St should be such that the following conditions [7] 


kySy 

< 1 

( 2 . 31 ) 

kzSz 

< 1 

( 2 . 32 ) 

uSt 

< 1 

( 2 . 33 ) 


are satisfied; where ky and kz are the wave numbers in the y and z directions respectively. If 
we assume ky = kz = ^ and for Sy = Sz = 6, the value of S must have to satisfy the following 
condition 



A 


6.28 


( 2 . 34 ) 


That is why we have selected ^ ^ for the FD-TD method. St is then chosen from the 

stability condition given by (2.22). 

An isotropic point source in two - dimensional case has a circular radiation pattern. We 
investigate how weU this pattern is maintained on the truncated mesh by using the Mur’s 
second order ABC. Figure 2.2 gives a clear view of the radiation pattern of the isotropic point 
source for two - dimensional case. It consists of the surface plot along with the contour plot 
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Figure 2.3: Radiatiou pattern of an isotropic point sovxce: mesli size 
no. of time steps = 650; surface plot and contour plot. 
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Figure 2.4: Contour plot of the radiation pattern after 650 time steps of an isotropic point 
source located at node (6,6) of a 76 by 76 node mesh. 

^also. The point source is located at node (36,36) and we present the result after 650 time 
steps. Though the number of time steps required to propagate the field to the boundary is 
only 80, we select 650 to show the effect of the reflection error introduced by this ABC. The 
error generated in one time step of the calculation accumulate and grow during the successive 
time steps. Prom the contour plot we see that the radiation pattern is almost circular with 
only sli^t deformations at the points near the boundaries and far away from the source. 
These errors are caused by the fact that waves with grazing incidence on a boundary axe 
not well absorbed, but partly (i.e. 10-15% depending on the distance from the source to the 
boundary) reflected. In Figure 2.4 we present the contour plot of the same dimensions except 
the source is located at the node (6,6). We see from the contour plot that now the patterns 
are not so circular as before and more deformations are introduced specially at the points far 
away from the source. That is why, it can be concluded that, to get more accurate results one 
should locate the source as far away from the boundary as possible. Since, (2.30) is very easy 
to implement and it gives exactly the same results as the more complicated one that follows 
from (2.27), we have used (2.30) for two - dimensional field problems. 

^If -we cut the siirface plot of a function of two variaHes by a plane paralld to the lower plane (i.e. horizontal 
plane) then the variation of the function on this plane is called the contour plot 

M 


Suppose, the point source is located at r* = {ys,Zs) = Now, this point source 

can be modeled by adding a term representing a current [5] 

hit) = (2.35) 

at Ts, where u(t) denotes the Heaviside unit-step function and C is a constant. That is, the 
input condition can be given by the following equation; 

= ^?”(io,A;o) + 4(i) (2.36) 

The above equation represents the isotropic point source at the node (jo, Aq) in the FD-TD 
method and C represents the amplitude of the electric field. 

2.4.2 Computed Results vs Theoretical Results 

Now, we will compare our computed results with the theoretical one. For the point source 
under consideration, the theoretical radiated field can be obtained by solving the equation: 

( d'^ I 3“^ \ d 

[d^ ~ " Mo^4(t)^(r-r,) (2.37) 

which is rotationaJly symmetric around and consists of the outgoing waves. The solution 
for the electric field distribution can be given as [10] 

E^ir,t) = -{^)l (i 

•X (2.38) 

The above equation can be given in the following simplified form: 

E^ir,t) = - cos - (i r - r-s p ^ /cq) 

X (lr-r,p+^2)-*d| (2.39) 

where | r - P . Equation ( 2.39) can not be put into a simple analytical 

formula. That is why, ( 2.39) is solved numerically using Simpson’s rule to get the theoretical 

value^ and the exact solution is plotted as a function of distance from the point source and 

^TMs is not the exact theoretical value because a very small amount of error is introduced by the Simpson’s 
rule. We have used Simpson’s 1/3 rule to solve the equation (2.39) and the error introduced by this method 

is given by Er = where m is the number of iteration and /*’' is the average fourth order derivative 

of the function which is integrated for the interval [0, ^i]. For our ca^ i = 650 6t and m = 800; and the 
maximum error introduced by the Simpson’s rule is in the order of 10 ^—10 This error is negligible and 
thererefore the results obtained by the Simpson’s rule can be considered as the theoretical value. 




Figure 2.5: Comparison of the theoretical results and computed results. Source is located at 
the node (10, 10). Results are compared for j — k nodes. 

at the same instant in time in Figure 2.5 . The results are compared for the diagonal nodes 
of the mesh(i.e. for j = k) that passes through the location of the point source. Point 
source is located at (10, 10) and no. of time steps = 650. Prom the Figure 2.4 we see that, 
computed results are very similar to the theoretical one and the error introduced is very less. 
A maximum 15 percent error is there which may be tolerated in some of the cases. 

2.5 Conclusion 

In this chapter we have given a brief review of the FD-TD method and the Mur’s second 
order ABC. We have studied the results using a point source only. We have also compared 
the FD-TD computed results for a point source problem with a more exact computanional 
result. In the next chapter we shall apply the FD-TD method for studying three coupling 
models as mentioned in the first chapter. 



Chapter 3 

Visualizing Optical Fields In 
Transmitter Packages 


In tMs chapter we present nninerical results about the optical field characteristics within 
the region between the source and the coupling end of the optical fiber for three coupling 
configurations. Since we are computing fields in a two dimensional domain, we approximate 
the fiber by a planar dielectric slab waveguide (Figure 3.1) for all the configurations under our 
consideration. This is composed of a dielectric slab of refractive index n\ ( =1.5 ) sandwiched 
between dielectric material of refractive index ( =1.499)^1, which we shall call a cladding. 
This represents the simplest form of the optical waveguide and can serve as a model to gain 
an understanding of wave propagation in optical fibers in a two dimensional domain. 

In the first configuration we consider a simple coupling model in which power is directly 
coupled into the waveguide — we call it as arrangement -1. In the second configuration we 
consider a diflferent type of coupling model in which power is coupled from the source into a 
taper-ended optical waveguide — we call it as arrangement-2. Finally, we consider a model in 
which power is coupled from the source through a microsphere into the optical waveguide — 
we call it as arrangement -3. Before studying the diflferent coupling configurations mentioned 
above we first study the propagation characteristics of a Gaussian beam in free space. We 
present results for two dimensional configurations using the FD— TD method and Mur’s second 
order ABC as before. A (top emitting semiconductor) laser is used as the optical source of 
wavelength A = 1.5)wm and in two dimensional configuration it is assumed as a line source. 
We present results for y-z plane; source is assumed along the y-axis at s = 0 and only the 
TE mode is considered as in chapter 2. For all the cases, the spot size of the laser beam 
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3.1 Diffraction Of A Gaussian Beam In Air 

Lasers emit Gaussian beams. In this section we present the study of the propagation of a 
Gaussian beam over a few micron distance away from the source in free space using the FD- 
TD method in two dimensional configuration. We assume only three field components; i.e. 
Ex{y,z,t), Hy{y,z,t) and Hz{y,z,t). At t = 0, the laser source of frequency / is assumed 
to be turned on. The diffraction of waves from this source is modeled by solving the finite- 
difference analog of the time-dependent Maxwell’s equations (i.e. 2.16, 2.20 and 2.21) of 
chapter -2. 

The software outline is almost the same as given in section (2.4) of chapter 2 except the 
source condition. In this case the laser is assumed to be a line source along y axis and located 
at A; = 10. Then the source condition for the FD-TD method can be given as follows [4]: 

10) = E^(j,10) + C exp sm(nu}St) (3.1) 

L ‘*’0 . 

where C is a constant and a; = 2 tt Co/A. At each point on the line^ the computer first 
calculates 10) in the normal manner of the algorithm and stores the value in memory. 
Then the value of the sinusoidal term is calculated and added to the stored value of E^{j^ 10). 
In effect, equation( 3.1) simulates the linear superposition of z directed plane wave and the 
ambient field along the grid line k = 10. Equation( 3.1) generates the desired sinusoidal 
incident plane wave. Note that, this source condition (i.e. eqn. ( 3.1) is also used for all 
models under our consideration. 

The field envelope, or the maximum absolute value, during the final half-cycle of time 
stepping is taken as the magnitude of phasor of the steady-state field field [4]. The intensity 
of the electric field ^ for the different values of z is plotted in Figure 3.2. From Figure 3.2 
we see that, with the increase in z^ the width of the beam is also increased; and it supports 
the theory of diffraction of a Gaussian beam in free space. In Figure 3.3, we have plotted 
the error plots for the above field intensity plots. In this figure, the computed results are 

^ First, tke absolute value of the electric field at each point in the computation domain is computed by the 
FD~TD method. This represents the mod of the electric field (i.e. Now, intensity (I) is given, by the 

relation I =; E#= |Ep. Therefore, the square of the envelope gives the intensity of the electnc field. 




y in meter 

Figure 3.2; Envelopes of the electric field for different values of compSO for z - 5A, complOO 
for ^ = lOA and comp200 for ^ = 20A; A = l.bixm . 


compared with the Fresnel diffraction formula [27]: 

j 2 Ctt 


Ex{y,z) = 


X(2z + jk oj^) 


exp < - j k 




2z (1 + w^u^lX^z^) 


oj{zy 


X exp 

where C is a constant and 

u{zy = 

and k = The intensity distribution can be given as 


2 (, , 


I{y,z) = A exp 


2y2 


u;izy 


where 


A = 


4 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


A^ (4z^ + fe^w^) 

which shows that the transverse intensity distribution remains Gaussian with the beam-width 
increasing with 2 : which implies diffraction divergence, wq is called the spot-size of the beam. 

This formula has less accuracy for small distaces from the source. That is why we see from 
the Figure 3.3, large amount of errors are introduced fox small distances from the source and 
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Figure 3.3: Error plots for the iiaussiau beam propagation in air. perrSO is for the electric fielc 
Z = 5A distance from the source, perrlOO for Z = lOA and perr200 for Z = 20A distai 
from the source, respectively. 
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Figure 3.4; Contour plot of the electric field for Gaussian beam propagation in free space. 

as the distance is increased error is reduced. In Figure 3.4, we have also given the contour 
plot of the electric field. From this figure we also see that, the width of the beam is increased 
with the increase in distance from the source and it agrees with the theoretical nature of the 
electric field. 

3.2 Study Of Model— 1 

In this section we study a simple coupling model in which optical power from the source is 
directly coupled into the fiber (Figure 3.5). Fiber is considered as a planner wave guide in 
y~z plane. The fiber under consideration is assumed to be a single mode step index fiber 
with core width = SAfim and cladding width = 60/^ as mentioned in the figure. The distance 
between the source and the coupling end of the waveguide is assumed as 45/im and the length 
of the waveguide is taken as 167/xm. 
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Figure 3.5: A simple coujiliiig model. Spot size of the Gaussian beam = nm. 


3.2.1 Software Implementation 

The main structure of the programming outline for model -1 is similar to that for studying 
the point source. In this case the dielectric constants for the nodes of the whole computation 
domain should be defined carefully. Some special treatments should be tahen for the nodes 
at the boundaries between two media. It can be proved that for the fields at the nodes on an 
interface between two media, the effective dielectric constant becomes for three media, 
it becomes gjid so on; where Ci, 62 , €3 etc. are the dielectric constants of the diflferent 

media. A proof of the above statement is given in Appendix A for a;n interface seperated 
by two media only. That is why, for the boundaries between the core-cladding, core-air and 
air-cladding, an average value of the permittivities is used. The source condition is the same 
as given by the eqn. (3.1). To apply the boundary conditions different phase velocity should 
be considered for different media. However, the programming outline for model— 1 can be 
given by the following steps. 


1. Begin: 

• INPUT PARAMETERS: 


rt/V 



' t . out ' 
0.675 
0.339 
0.00192 
-0.335 
-0.671 



0 50 100 150 200 250 


Figure 3.6: Contour plot of the electric field between source and fiber-end for model -1. 





2 . Read: miinber of nodes in air on both sides of the waveguiec in y-direction {yomax), 
number of nodes in cladding on both sides of the core iyimax), number of nodes in 
the core region in the y-direction {y 2 max)-, number of nodes between the source and 
the coupling end of the waveguide in z-direction (zimax), no. of nodes in fiber in z- 
direction [z-zmax), operating wavelength (A) of the source, dielectric constant of the 
air (eo), dielectric constant of the core (ci), dielectric constant of the cladding (€2) and 
permeability of the media (/io)- 

• SELECTION OF SOME PARAMETERS FOR THE FD-TD METHOD: 

3 . Calculate: phase velocity in the air(upo) = pbase velocity in the core (upi) = 

and phsLse velocity of the cladding {Vp2) = 

4 . Select: maximum phase velocity (vp^ax) = ^po = Wia® < Vpi,Vp2,Vpo > 

5 . Select: Sy = 6 z = 6 = A and St = 5— ^ — 

6. Select different dielectric constants for the different nodes depending on the region in 
which the nodes are located. 

7 . Take an average value of the dielectric constants for the nodes at the interfaces between 
different media. 

8. Select: no. of time steps (nma®)=inaximum of the no. of time steps required to reach 
the wave at the boundary in the y-direction from source and no. of time steps required 
to reach the wave at the boundary in the z-direction from the source. 

• STARTING OF THE CALCULATION OF THE FIELD COMPONENTS: 

9 . Start time stepping (n) at n = 0 

10. Choose a line along y-axis (i,^o) (say). Apply input along this line as given by the 
equation 3.1. Source should be located as far away from the boundary as possible. 

• CALCULATION OF THE FIELD COMPONENTS AT THE INTERIOR NODES: 

11. Calculate Hy(j, k), H^{j, k) and k) at the interior nodes of the computation domain 

with the regular finite difference scheme. 

• SPECIAL TREATMENTS FOR THE BOUNDARY NODES: 



12. Apply Mur’s second order ABC to calculate the electric field components at the bound- 
ary nodes. Different phase velocities are considered at different nodes on the boundaries 
for different media. 

13. Increase n by one i.e. n — *■ n -f 1 

14. Verify: if n < n,„ax then go to step ( 9); else finish, 
o CALCULATION OF THE FIELD ENVELOPE: 

15. Compute the time period (T) of the sinusoidal input. Calculate the no. of time 
steps (Uenv) half cycle (T/2). 

16. Find the maximum of the field components at each node during the final Uenv time 
steps. This will give the envelopes of the corresponding field components. 

17. Print output. 

18. End 

3.2.2 Results 

First, we present the result for the field characteristics between the source ajid the coupling 
end of the waveguide. In Figure 3.6 we present the contour plot of the electric field within 
the region between the source and the fiber-end considering the model given in Figure 3.5. 
The contour plot for the electric field is already given when there is no fiber i.e. propagation 
of Gaussian beam in free space, in Figure 3.4 for the same computation domain as between 
the source and fiber-end. If we compare the two contour plots we can see the differences 
easily. From Figure 3.4 we see that the phase fronts are very regular when the Gaussian 
beam propagates in free space. But from the Figure 3.6 we see that phase phase fronts are 
not regular and a weak interference pattern is generated. This is mainly due to the reflection 
ocurred at the air-core and air-cladding interfaces. 

In Figure 3.8 we present the electric field profiles taken at three positions within the source- 
fiber region. Prom these profile plots we see that the difference between the field profile at 
Z — 7.5/im and Z = 18/im is very less but there is a significant difference between the above 
two and the plot at Z = ZOfjtm. These are the envelope plots which show the variation in 
amplitudes only. 
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Figure 3.8: Field profiles within the source-fiber region for naodd-l. usa.ll for 2, - ISpim; 
usa.l2 for Z = ISfini’, and usa.l3 for Z = SOfiTfi from the source, respectively. 



Now we present the fundamental field profile within the waveguide and the coupling effi- 
ciency. A theoretical treatment of the fundamental mode profile for TE mode in a symmetric 
planar waveguide is given in appendix B. In Figure 3.7 we present the fundamental field pro- 
files within the waveguide at the distances Z = 175/i and 205nm from the source along with 
the theoretical value. From the figure we see that the field profile within the core is similar to 
that of the theoretical one. But within cladding region the field profile deviates to some large 
extent from the theoretical value. This is mainly due to the short length of waveguide under 
consideration. Over which the evanescent fields in the cladding region have not settled down 
to the final equilibrium values. Usually it takes about 50m distance for the evanescent fields to 
the equilibrium value resulting resulting in propagating modes. It is practically impossible to 
simulate over 50m long by the FD-TD method. For the field profile obtained at Z = 205//m 
is used to calculate the coupling efficiency of this model. The coupling efficiency obtained 
by this model is only 38.5027%. The maximum obtainable theoretical coupling efficiency is 
65.23% (see Appendix C). This dissimilarity is mainly due the planar structure consideration 
of the fiber. Moreover, we have considered a very short length of the fiber which is practically 
impossible to set up a particular mode within the fiber. Further, for the theoretical consider- 
ation cladding is assumed to be extended to infinity whereas a finite value is considered for 
our study. A small error is also introduced by the FD-TD method. 

3.3 Study Of model-2 

In this section we study a different type of coupling model in which taper-ended fiber is used 
to couple power from the source. One of the models of this type is shown in Figure 3.9. The 
dimensions and the refractive indices are mentioned in the figure. 

3.3.1 Software Outline: An Extra Treatment 

The main programming outline for model-2 is exactly the same as for model-1 except the 
tapper portion of the fiber. The extra care is taken only to define the dielectric constants 
for different nodes within different regions carefully. In this case the extra tapper portion is 
discretized using stair case approximation. The input parameters and other things are same 
as for model-1. That is why in this section we present only the discreatization procedure 
of the tapper portion of the fiber. The tapper portion is shown in Figure 3.10. To define 
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Figure 3.9: Model-2: coupling witb taper-ended fiber. Spot size of tbe Gaussian beam —nm. 
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Figure 3.10: Tapper portion of the fiber of model-l. 


dielectric constants for the region ADEB is very easy. But, for the regions ABC and DEF, we 
use the following method. Suppose P {y,z) is a point located anywhere in the left-side region 
of CF. The slope of the line AP is given as 


TO = tan(^) = 


z - Zl 

y-yi 


Now the slope of the line AC is fixed for a particular structure. This slope is given as 


TO/ = tan(^i) = — 

^ y2-yi 

If TO > TO/, the point P is located within the air; otherwise it is located within the triangle 
ABC(cladding). The nodes at the interfaces are selected in the following. From the point P 
draw the perpendicular PN on the line AC. Let the angle between the lines AC and AP be 
Od- Now, if a point satisfies the constraint 

PN = AP siniOd) < 2 


we consider this point as the interface point between the air and the core for the stair case 
approxinaation of the interface; and we assign an average dielectric contant for this point. 



UsiEg the above idea we can select the nodes within the core, cladding and air. The same 
rule is also applied for the region DEF. After selecting the nodes for different regions we can 
eissign different dielectric constants for different nodes. For the nodes at the interface between 
different media we use an average dielectric constants as mentioned earlier. 

3.3.2 Results 

Now we present the results for model-2. In Figure 3.11 we present the contour plot of the 
electric field within the region ABCD as shown in the figure. From the contour plot we see 
that a strong reflection occur at the air-core and air-cladding boundaries. As the angle of 
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incidence of wave is different for the different points of the taper portion of the fiber, different 
amount of energy is reflected from different sections. Most of the energy is reflected from the 
taper portion of the fiber. Due to the normal incidence of the field on the core re^on, a small 
portion(maximum 4%) of the energy is reflected. Therefore, a narrow field is propagated 




through the core of the fiber. Source condition and air-core, air-cladding and core-cladding 
boundaries are treated in the same way as before. In Figure 3.12 the fundamental field profile 
for Z = 205fj,m from the source is plotted along with the theoretical plot. If we compare 
Figure 3.7 and Figure 3.12 we see that field profile for model-2 is more accurate both in 
core and cladding to the theoretical value than model-1. The coupling efldciency obtained 
by this model is 51.937%. In this case we see that there is an improvement of the coupling 
efficiency over the butt coupling. This is due to the fact that a narrow field is entered into 
the fiber and the taper portion acts somehow as a lens. In practical cases coupling from 
semiconductor laser into single mode fibers is done with tapered hemispherical ends[21]. The 
taper-end of the fiber collects most of the energy in the core region and that is why a strong 
field is propagated through the core and increase the coupling efiiciency. 

3.4 Study Of Model— 3 

In this section we study a good model using a nonimaging microsphere[23,24] between the 
optical source and fiber as shown in Figure 3.13. First we make some assumptions: the sphere 
lens has refractive index ni = 2.0 and outside medium is air (refractive index riom = 1-0). 
Then the focal length can be calculated from the gaussian lens formula [30] 

^ __ ^om /g g\ 

s q T 

where s and q are the source and image distances, respectively, as measured from the lens 
surface, and r is the radius of curvature of the lens surface. The following sign conventions 
are used to solve the eq. ( 3.6): 

• Light travels from left to right. 

• Object distances are measured as positive to the of a vertex and negative to the right. 

• Image distances are measured as positive to the right of a vertex and negative to the left. 

• All convex surfaces encountered by the light have a positive radius of curvature, and concave 
surfaces have anegative radius. 

Using these conventions we now find the focal' point for the right-hand surface of the lens 
shown in the Figure 3.13. To find the focal point, we set q = oo (for collimated output) and 
solve for s in eqn. 3.6, where s is measured from point B. With r = —Df2, eqn. 3.6 ymlds 
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Figure 3.13: Nonimaging microspliere is placed between the source and the fiber; model-3 
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S = f = D. Thus the focal point is located on the lens surface at point A. This, of cource, 
changes if the refractive index of the sphere is not equal to 2.0. 

That is why for n; = 2.0 the source should be placed close to the surface of the lens and 
fiber can be placed right side of the lens surface. Generally the fiber is also placed as close to 
the lens as possible for the compact packaging purpose. 

The use of homogeneous glass microsphere lenses has been tested in a series of several 
hundred laser diode assemblies by Khoe and Kuyt [26]. Spherical glass lenses with a refractive 
index of 1.9 and diameters ranging between 50 and 60/im were epoxied to the ends of 50/^m 
core-diameter graded-index fibers. Coupling efficiencies in these experiments ranged between 
50 and 80 percent. 

3.4.1 Software Outline: An Extra Treatment 

In this case also the main programming structure is the same as model-1 except the lens 
between the source and waveguide. Extra treatment is required only to define the dielectric 
constants for the nodes in the lens material instead of the air between the source and the 
wavegmde. In two dimensional configuration the lens is assumed to be a cylindrical with a 
circular cross section on the y-z plane. Therefore we discuss the extra treatments required 
for the lens. To select a node inside or outside of the lens material (i.e. in air region) we use 
the rule as follows. The equation of a circle in y-z plane can be given as 

y^-\-z^ = r^ 

where r is the radius of the circle. Let P(y,z) be any point within the region between the 
source and the waveguide. Now, if 

y'^ + z^ -r^ <0 

then the point is located within the circular region; otherwise it in the outside of it (i.e. in 
air). For points at the interface between the lens material and the air we use the following 
condition 



Now we can select the location of the different nodes and assign different didectric constants 
depending on the medium in which it is located. For the interface nodes an averaging rule of 
the dielectric constants is used. 





3.4.2 Results 


We have studied the lens with refractive index m = 2.0, diameter D = A2{im ;and for a 
step index fiber by the FD-TD method. Fiber is placed just 1.5/im right of the point A. 
The input condition is the same as given by 3.1. In Figure 3.14 we present the contour 
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Figure 3.14: Contour plot of the electric field for model-3 

plot of the electric field to see the field pattern generated between the source and the fiber. 
Prom Figure 3.14 we see that as the field propagates towards the fiber coupling end a part 
of the input energy is reflected at the lens-air interface. A part of the energy is also reflected 
from the air-core and air-cladding interfaces. Maximum 30% of the input field amplitude is 
reflected for normal incidence but different amount of energy is reflected from different points 
of the lens surface due to different angle of incidence at different points. The most important 
thing is to note that inspite of the reflection occured at the interface of the air and the lens 
material, a strong field is propagated through the lens material towards the fiber along the z 
axis and this actuUy shows the improvement of the coupling efiicieiicy of this lensing scheme. 





The fundamental field profile is presented in Figure 3.15. From this figure we see that again 
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Figure 3.15: Fundamental field profile at Z = 205/im from the source for model-3. 

the field is not stable in the cladding region but the field within the core is almost stable. We 
also observe from Figure 3.15 that in this case most of the energy is confined within the core 
region. For this model the efficiency obtained by the FD-TD method is 69.394%. In this case 
we see that the coupling efficiency is the highest among all the three coupling attachments. 
This is mainly due to the the magnification of the emitting surface area of the source by 
the lens [29, chapter-5]. The field profiles in between the lens and the fiber is shown in the 
Figure 3.16. We see from the figure that a narrow field enters into the fiber. If there is no 
lens, due to the diffraction effect the width of the light beam spreads out resulting in a small 
amount of energy incident on the core of the fiber. But, owing to the presence of the lens 
light beam does not spread out and therefore a strong field is propagated through the fiber. 
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Figure 3.16: Field profiles in between tbe lens and fiber for niodel-3. 12us.ol: for no sWfting; 
12us.ll for shifting the source 3.0 fj/m in the downward direction; and 12us.21 for shifting the 
source 1.5 um in the upward direction. 




That is why in this case the coupling efficiency is much higher than the other two cases. 
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Figure 3.17: Contour plot of the electric field for 1.5^m shifting of the source in the downward 
direction laterally. 

Now we present the results for the misalignment problems. From the field profiles in 
between the lens and the fiber end face presented in Figure 3.16 for the shifting of the source 
in the downward and upward directions, we see that due to the shifting field strength is 
increased before the cladding regions. Moreover, the light beam spreads out due to the 
shifting. That is why, less energy is propagated through the core of the fiber. The contour 
plot of the electric field within the source-fiber region is presented in Figure 3.17 when the 
source is shifted 1.5/um in the downward direction laterally, if we compare Figure 3.14 with 
Figure 3.17 we can easily see that less amount of energy is propagated through the lens 
towards the fiber. In Figure 3.18 we present three field profiles ai Z = 205/J’m from the 
source: 12us.o3 for no shifting, 12us.l3 for the source shifted 3pm in the downward direction 
and 12us.l3 for the source shifted l.hpm in the upward direction. Due to the shifting of 
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Figure 3.18: Field profile plots for tlie misalignment problem for model-3. 12us.03: for no 
shifting of the source; 12us.l3: for shifting the source 3.0 jum in downward direction; and 
12us.23: for shifting the source 1.5 in the upward direction. 



the source, we see that the amplitude of the field is reduced significantly. For 12usl3 and 
12usal3 the coupling efficiencies are 33.436% and 54.193%, respectively, while for no shifting 
this is 69.394%. Therefore, the coupling efficiency is reduced for the shifting of the source 
significantly. Further the power propagated through the cladding region is increased in this 
case and cladding fields are not in equilibrium state due to the short length of the fiber under 
consideration. 

3.5 Conclusions 

In this chapter we have given the field characteristics between the source and the planar op- 
tical waveguide for three coupling attachments. We have also studied the coupling efficiencies 
for different coupling models. We see that for model-1, the coupling efficiency is 38.5027% 
which is large enough compared to the practical coupling results for actual fibers. For actual 
fibers the coupling efficiency lies between 10 and 20 %. For model-2 we get the coupling 
efficiency 51.937% which supports the practical results. Practically, the power coupling from 
a semiconductor laser into single mode fibers with hemispherical ends gives the coupling ef- 
ficiency between 40% and 70%. We see from model-3 that coupling efficiency is the highest 
among all the three coupling attachments under our considerations. In this case we get the 
coupling efficiency is 69.394%. For the actual fiber this efficiency lies between 50% and 70%. 

The first configuration is the simplest one. It is very easy to handle for practical experi- 
ments. The second coupling arrangement is sensitive to the proper mechanical alignment and 
creat additional complexities for the practical cases. The third coupling arrangement creats 
more complexities than the other two. One problem is that the lens size is similar to the 
source and fiber core dimensions, which introduces fabrication and handling difficulties. But, 
for the efficient coupling the last coupling arrangement is the best among the three as it gives 
the highest coupling efficiency. 

We have also developed the visualisation soft-ware package to see the optical fidd char- 
acteristics within the source-fiber region of the transmitter packages. The advantages of this 
package are: (1) we can change the parameters i.e. the refractive indices of the core, cladding 
and the medium within the source-fiber region, length of the fiber, wavelength of the source, 
distance between the source and fiber, dimensions of the core and cladding regions, source 



emitting area etc. can be changed as these are given as the inputs to the programming; (2) 
The absorbing boundary condition is very simple and easy to implement. (3) Soft-ware is 
written in FORTRAN 77 which is also easy to handle and understand. Using this soft-ware 
we can study the reflection characteristics of the field from the fiber end and the coupling 
eflSciency for different transmitter packages. Varying the input parameters it is possible to 
find out a particular arrangement condition for a particular transmitter package for which an 
efficient coupling can be achieved. Following this optimum coupling condition, one can de- 
sign this particular transmitter package for the practical purpose so that maximum coupling 
of power can be obtained. Moreover, this soft-ware can be used to study the effect of the 
mechanical misalignment for the fiber optic transmitter packages. 



Chapter 4 

Further Scope Of Work 


la tMs chapter we present the further scope of work for which our results and discussions may 
be used as the initial required background. Using the FD-TD method and Mur’s ABC, we 
have studied the optical field characteristics within the region between the source and the fiber 
for only three coupling configurations. Following the same approach, we can study the field 
charateristics and the coupling efficiencies for other coupling configurations as suggested in the 
first chapter. Moreover, the results can also be extended for three dimensional configurations. 
But, in this case the major difficulties are the requirement of a large amount of computer 
memory and computation time. One may study the coupling problem for the three main 
types of lenses discussed in the following section. 

4.1 Three Main Lenses For Coupling 

Now we present three main types of single-element lenses for fiber-optic devices; (1) quarter- 
pitch GRIN-rod lens; (2) ball lens and (3) rod lens. These three lenses are very useful 
fiber-optic devices which can be used for the optical power coupling from the source into the 
fiber and from fiber to fiber. 

4.1.1 GRIN-Rod Lens: 

The refractive index in a GRIN-rod lens is given by 

n'^{r) = n§[l - (g rf + {g rf -h he (g r)®], 

where n(r) =refractive index at distance r from the axis; 
no =Tefractive index on the axis; 
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Figure 4.1: Fundamental imaging parameters of the GRIN-rod lens. 

g =focusing constant; and 
hiyhe =higher order constants. 


For paraxial rays the higher order terms with coetRcients can be neglected. A ray incident 
on a GRIN-index lens travels along sinusoidal curve with a periodic length (pitch)= y. For 
the usual quarter-pitch lens, the focal length and the position of the principal points H are 
(Figure 4.1) 
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4.1.2 Ball Lens: 

A ball lens is fully described by its refractive index and radius r. Application of the thick 
lens formula leads to (Figure 4.2) 

H = r. 


4.1.3 Rod Lens: 


One focal point of the rod lens is chosen at the plane face of the lens, so that only independent 
paxameters are refractive index n and radius of curvature r (Figure 4.3). Using the thick lens 


central L'8RARY 

I I.T., <ANPUR 





formula we get. 


4.1.4 Complete Lens System: 



For a maximum coupling between the source and the fiber or between the fibers, an one to 
one image should be made of the emitting surface upon the receiving one. For fiber to fiber 
coupling, special consideration should be taken so that the launching angle with the fiber end 
face has to be the same at the receiving fiber end. All these conditions axe satisfied only if 
the ends are positioned in the focal plane while the foci of the lenses coincide (Figure 4.4). 
The incoming light cone is tilted if the foci of the lenses do not coincide. In a proper design, 
this tilt angle should always be small compared with the numerical aperture of the fiber. 

In a well designed lens coupling each currently proposed lens type i.e. GRIN-rod lens, 
ball lens, and rod lens, will reach equal imaging quality and hence equal couple efficiency. For 
typical communication fibers, coupling efficiencies up to 95% axe possible. 



4.2 Study With Antireflection Coatings: 

In this section we present the idea of studying with antireflection coatings applied to optical 
surfaces in order to modify their reflectance/transmittance properties. About 4 % of visible 
light is reflected at any uncoated air-to-glass surface of a lens, prism, or other optical com- 
ponent. For designs using more than a few components, losses in light level can accumulate 
rapidly. More important are the corresponding losses in image contrast or modulation due 
to weak reflected ’’ghost” images superimposed on the desired image. Such unwanted images 
are often defocused beyond recognition so that contrast reduction (rather than confusion) is 
their primary effect. 

We are able to apply coatings designed to operate within a wavelength region from 200 
to 1600 nm. Coating types include single layer dielectrics, multilayer dielectrics and hybrids. 
Some coatings are for use at a specific wavelength while others cover a broad spectral range. 

4.2.1 Single Layer Coatings: 

Surface reflectance depends on the difference in refractive indices for air and glass. For visible 
light this is approximately the difference between 1.0 and 1.5. By selecting a material with an 
index of refraction that falls between these two values (near the geometric mean), and applying 
a coating of approximately 1/4 -wavelength optical thickness to the surface, it is possible to 
reduce the reflectance significantly. Magnesium fluoride (M^Fz), a durable substance with 
a thin film refractive index of 1.38 at 550 nanometers, is the most commonly used coating 
material for single-layer antireflection (AR) coatings. 

Instead of one reflection at the air-glass interface, there is a first reflection at the air- 
M 5 F 2 interface followed by a second reflection at the MjF 2 -glass interface. These reflections 
interfere destructively, resulting in a minimum reflectance within the visible spectrum at 
normal incidence of 1.5% or less, as compared with with the original 4% for the uncoated 
glass surface. From the principle of energy conservation, confirmed experimentally, there 
must be a corresponding increase in overall transmission of a coated surface, including both 
interfaces. 

Such coatings work extremely well over a wide range of wavelengths and angles of inci- 
dence., despite the fact that the theoretical idea of zero reflectance is achieved by a film of 



1/4- wavelength optical thickness only for normal incidence, and only if the refractive index 
of the coating material is exactly the geometric mean. In fact, the single layer 1/4- wave 
thickness MgF^ coating designed for normal incidence makes its most significant contribution 
to the transmission of steep surfaces, where most rays are incident at very large angles. 

4.2.2 Multilayer Coatings: 

By combining multi-layers of high, medium and low index materials it is possible to cause 
the reflected and transmitted beams to combine selectively in order to achieve an important 
performance over a single layer coating. The improvement can be in broadening the spectral 
bandwidth of the coating or in almost nullifying the reflectance at a particular wavelength. 

The typical High Efiiciency BroadBand AntiReflection (HEBBAR) coatings give a very 
low reflectance over a broad spectral bandwidth. The normal (0°) coating typically gives less 
than 0.5% reflectance over a range of 440 to 660 nm. This coatings can be modified for use 
at 45® or can be shifted to cover a wavelength range from 750 to 1100 nm. 

To obtain the ultimate in transmittance at a single wavelength, it is possible to deposit 
a ”V-coat” which will reduce the reflectance to less than 0.25%, and typically closer to 0.1% 
for that one wavelength at normal incidence only. 

4.3 Mismatching Study: 

When power is coupled from the source into the fiber or from fiber to fiber, the center of the 
emitting area of the source or of the transmitting fiber should be located on the aods along 
the receiving fiber length to get the maximum coupling efiiciency. We have studied the effect 
of mismatching between the source and the receiving fiber by shifting the source laterally; 
and calculated the coupling efficiencies for only two specific cases. One can extend this idea 
to study the effect of mismatching on the coupling efficiency by giving lateral and angular 
displacements to the source or the emitting fiber for both the source-fiber and fiber-fiber 
coupling problems. This study may be very useful for designing the transmitter packages, 
splice and connectors, respectively. 



4.4 Remarks: 


In this chapter we have presented only some elementary ideas about the further scope of 
works. These studies may be very useful for the optimum design of the optical systems. 
These ideas can also be applied to the different type of fibers with different kind of refractive 
index distributions. 



Chapter 5 

Conclusions 


The main purpose of this thesis is to visualize the optical field characteristics in transmitter 
packages. In this regard we have studied three different coupling attachment schemes named 
model— 1, model— 2 and model— 3. We have also studied the the coupling efficiencies for these 
three coupling arrangements. Finally, we have studied the misalignment problem for model-3 
only. In this case the source is shifted laterally in the upward and down-ward directions; and 
we have studied the effect of the coupling efficiency and the field patterns for this type of 
simple misalignment problem. The whole results are given as follows. 

• We have studied the FD-TD method and Mur’s second order absorbing boundary con- 
dition. We have presented the results for a point source problem to find the accuracy of the 
FD-TD method. The computed result is compared with the theoretical result (which is also 
obtained by a more accurate numerical method; i.e Simpson’s rule; for the point source prob- 
lem), and we have got an error of maximum 15%. Average dielectric constant is taken for the 
boundary nodes and Sy = 6z = A/10 is considered to implement the FD-TD method. The 
FD-TD method is implemented only in two dimensional configurations (i.e. in y-z plane). 

• In model-1 we have studied the optical field patterns within the source-fiber region and the 
coupling efficiency for a flat fiber end face centered directly over the source. We have pre- 
sented the" reflectio n cEajacteristics of the optical field from the fiber end face by the contour 
plot. The fundamental field profiles are also given. The coupling efficiency obtained in this 
case is 38.5027%. 

• A taper-end fiber is used in model— 2 instead of the flat fiber end face; In this case we have 
presented the field characteristics by a contour plot. The propagation of the dectric fidd from 
the source towards the fiber is shown in this contour plot. The reflection characteristics of the 



field from the fiber end is also shown in this plot. We have presented the fundamental field 
profiles within the fiber for this model. Finally, we have calculated the coupling efficiency. 
The coupling efficiency for this model is 51.937%. 

• An efficient coupling arrangement is studied in model— 3 by placing a microsphere lens in 
between the source and the fiber. In this case we have presented the reflection of the electric 
field from the lens material and the fiber by a contour plot. Prom this plot we have seen that 
astrongfield is propagated through the lens and propagates into the fiber. The fundamental 
field profile within the fiber is also shown for this model. In this case the coupling efficiency 
is 69.394%. Finally, we have studied the misalignment problem for this model by shifting 
the source the source laterally in the upward and downward directions. The effect on the 
reflection characteristics of the field is shown by another contour plot. It is also seen that due 
to this misalignment the coupling efficiency is reduced to a significant amount. 

We have studied these problems by the FD-TD method in two dimensional configurations. 
The fiber is assumed as a slab waveguide. We have developed the required soft- ware in 
FORTRAN 77 computer programming language. The main difficulties for the FD-TD method 
is the requirement of large amount of computer memory and computation time. To study 
each model the approximate execution time is 25 hours in DEC-ALPHA machine at IIT 
Kanpur. 
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Approximate and can be given as 


dz Ji 


Hyjm) - Hy{m - 

Sz 

2 

Hy(m + i) 


where m is assumed to be the position of the interface, and m + ^ and m- \ denote the 
positions a half step above and below the interface, respectively. 

Substituting ( A. 6) and ( A. 7 ) into ( A. 4 ) and ( A. 5 ), we get 

Hy{m) « - ] Hy{m + -) 

€1+62 z 

, ^2 tt / 1 62 - €1 dHz 6z 

+ ^ — fA.8) 

61+62 2^^ 61 + 62 % 2 ^ ’ 

Substituting the Hy^m) value of ( A.8) back into ( A.6) and ( A. 7 ), substituting the 
resulting ( A.6) and ( A. 7 ) into ( A. 4 ) and ( A. 5 ), and adding together, we get 

61 + 62 dEx _ dHz _ Sy{m + I) - Hy(m - |) 

2 dt ~ dy Sz 

- fA9) 

“ dy Az 


Comparing ( A.2) and ( A.9) we can conclude that for the points at the interface between 
two media, the effective dielectric constant should be taken as the average of the dielectrics 
of the media. Similarly, it can be proved that, this averaging rule is also applicable for more 
than two media. 



Appendix B 

TE Modes Of A Symmetric Step 
Index Planar Waveguide 


Now we carry out a brief modal analysis of a symmetric step index planar; waveguide. We 
consider only TE mode with Ex, Hy and components ; other components are assumed to 
be zero. If we choose the z-axis along the direction of propagation, then without any loss of 
generality, we can write the components in the form of 


Ex = Eoxiy) exp[j(wt - /3z)] 

(B.1) 

Hy = Hoyiy) exp[j(wt - ^z)] 

(B.2) 

Sz = Hoz{y) exp[j(wt - ^z)] 

(B.3) 

where we have assumed that field does vary in the x-direction; Eox{y), Soy{y) and Hoziy) are 

the arbitrary field profiles which depend only on y, and ^ is the propagation constant. Now 

from Maxwell’s equations we have 


Hy = -^Ex 

(B.4) 

TT _ i 

ixou dy 

(B.6) 

je,n\y)Ex = j0Hy + ^ 

(B.6) 


where n^{y) is refractive index distribution which is assumed to be independent of a: and z. 
Substituting ( B.4) and ( B.5) in ( B.6) we get the following key equation 


^ + lkln\y)-f]Ex = 0 (B.7) 
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where 


ko = wy^eo Mo 
u 



R.I.= rii 

R.I.= 712 


Cover (Cladding) 



Figure B.l: (a) The simplest planar/ optical waveguide (b) The refractive index distribution 
for a symmetric fpfemar waveguide 

Our analysis is valid for an arbitrary y dependent profile. We now assume a specific field 



profile given by (see Figure B.l) 


n(y) = I "'i W < '*/2 

1 « 2 ; \y\ > rf/2 


(B.9) 


where Ui and n 2 are the refractive indeces of the core and cladding, respectively and d is 
the width of the core. To solve the eqn. ( B.7) the boundary conditions must be satisfied 
at the discontinueties. The boundary conditions can be stated as: and ^ must be 

continueous at a: = i:d/2. Substituting for n{y) in eqn. ( B.7) we get 


+ (^0 " Pj Ex = 0; \y\ < d/2 film 


(B.IO) 




fiy 2 + (^0 «2 - r) = 0; |j/| > d/2 cover (B.ll) 

For guided modes it requires that the field should decay in the cover (in the region y > d/2) 
and thus we must have 


> ^0 *^2 

Furthermore, to satisfy the boundary conditions at y = ±d/2, we must have also 

< kl n\ 


(B.12) 


(B.13) 


That is why for guided modes we must have 


ko nl < < k'o n( 


,2 ^2 


(B.14) 


Therefore we can write the eqns.( B.IO) and ( B.ll) in the form 


where 


and 


dy"^ 

dy^ 


+ = 0; |j/l< d/2 film 

(B.15) 

+ 'y'^Ey = 0; |y|> d/2 cover 

(B.16) 

= klnf- 

(B.17) 

7^ - p^~klnl 

(B.18) 





The solution of eqn.( B.15) can be given in the form 


Ex(y) = A cos{y) + B sin(K y); for \y\ < d/2 (B.19) 

where A and B are constants. The solution of the eqn. ( B.16) can he given in the form 

„ f . f Ce-^y; y< -d/2 

= [oe-^y, y>d/2 

where C and D are also some constants. A, B, C and D are determined from the boundary 
conditions. Now if we assume that the refractive index distribution is symmtric about y = 0, 


i.e. 


T?{~x) = n^(x) 


(B.21) 


then the solutions are either symmetric or antisymmetricfunctions of y\ thus we must have 


Exiy) = Ex{y) symmetricmodes 
Exiy) = Ex{y) antisymmetricmodes 

For the symmetric mode we must have 

P _ I AcosiKy)] \y\ < d/2 
• - \ Ce-'y Ivl; \y\ > d/2 

Using the continuity of Ex{y) and ^ at » = ±d/2 we have 

A cos(k d/2) = C e”'’’ 


(B,22) 

(B.23) 


(B.24) 


(B.25) 


and 


— K A sin(K d/2) = C e 


(B.26) 


respectively. Prom ( B.25) and ( B.26) we get 


Since 


(«d/2) tan(Kd/2) = d/2) 

^2 _ ^2 _ ^2 ^2 
= feo(”i - ~ 


(B.27) 


(B.28) 



we can write 



COtf 



(B.34) 


Eqns. ( B.32) and ( B.34) are very useful equations to calculate tlie propagation constant 
/? of the modes. First, for a particular value of V ( B.32 for symmetric mode) or ( B.34 for 
antisymmetric mode) is calculated numerically to obtain the value of Then is calculated 
using ( B.30). For example, if ni = 1.5, 712 = 1.499 and d = 8.4/im from eqn. ( B.31) 
we get V = 1.927661. Then using eqn. ( B.32) and the Newton- Raphson as the numerical 
method we get ^ = 0. 722805 (for symmetric mode). Finally, from eqn. ( B.30) we calculate 
/? = 6.283357 X 10®. 7 = 0.1518075 x 10® and k = 1.720964 x 10® are also calculed from 
eqns. ( B.29) and ( B.30), respectively. 

Now we give two important points about the modes. 

• If0<V <7 rwe have only one discrete (TE) mode of the waveguide and this mode is 
symmetric in y. For this case the waveguide is called to be a single moded waveguide. For 
example, if A = l,5)tim, ni = 1.5 and n2 = 1.499 then for a single mode operation we must 
have (from eqn. (B.31) 




. Solving we get 


d < 13.3)LtTO 


• In general, if 


2mw < V < (2 m+ l)7r 


there will be (m + 1) symmetric modes, and m antisymmetric modes and if 


(2 TO + l)7r < F < (2 m + 2)ir 

then there will be (to + 1) symmetric modes, and (to + 1) antisymmetric modes where to = 

0) 1? 2, 3, The total number of modes will be the integer closest to (and greater than) V/tt. 

Therefore, for ni = 1.50, = 1.499, A = l.SO/im and d = 8.4/iTO, we have V = 1.927661 and 

therefore only one mode (symmetric) exists within the waveguide. 



Appendix C 


Coupling Efficiency For Butt 
Coupling 


For single mode fiber the fundamental field profile can be approximated to a Gaussian beam 
with spot size wj which is given by 


Wf 


d 0.65 "I" 


1.619 2.879\ 

71.5 + J 


(C.l) 


where F is a dimension less parameter of the fiber (see Appendix B) and a is the radius of 
the core. For butt coupling i.e. centering the fiber end face directly over the (laser diode) 
source, the maximum coupling efficiency can be given by 


V = 



4 



ko[wl+w^j) 


(C. 2 ) 


where Wq is the spot size of the laser beam, z is the distance between the fiber and the source 
and ko = 27r/A. 

For model-1 V = 1.92 (see Appendix B), Wo = 3.5fJ,m and z = 45jum. From equation (C.l) 
we have Wf = 5.62727/tm. Now using equation (C. 2 ) we can calculate the coupling efficiency 
for A = 1.6/im as 77 = 0.6523 = 65.23%. 

Generally, for the butt coupling case coupling efficiency ( 77 ) is very less (less than 15%). 
As Wf approaches to wq the coupling efficiency is approaches to 100% for z = 0. For the 
practical cases Wo < 1.0. But for our model - 1 we have considered wq = 3.5 fim which is very 
high and that is why the coupling efficiency is also very large compared to that obtained by 
general experiments. 
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